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When vortices are displaced in Bose-Einstein condensates (BEC), the Magnus force gives the 
system a momentum transverse in the direction to the displacement. We show that Bose-Einstein 
condensates (BEC) in long channels with vortices exhibit a quantization of the current response 
with respect to the spatial vortex distribution. The quantization originates from the well-known 
topological property of the phase around a vortex it is an integer multiple of 27 t. In a similar way 
to the integer quantum Hall effect, the current along the channel is related to this topological phase, 
and can be extracted from two experimentally measurable quantities: the total momentum of the 
BEC and the spatial distribution. The quantization is in units of m/2h, where m is the mass of 
the atoms and h is Planck’s constant. We derive an exact vortex momentum-displacement relation 
for BECs in long channels under general circumstances. Our results presents the possibility that 
the configuration described here can be used as a novel way of measuring the mass of the atoms in 
the BEC using a topological invariant of the system. If an accurate determination of the plateaus 
are experimentally possible, this gives the possibility of a topological quantum mass standard and 
precise determination of the fine structure constant. 

PACS numbers: 03.75.Lm,73.43.-f,06.20.fb 


I. INTRODUCTION 


One of the striking aspects of the integer QHE is the 
very precise quantization of the transverse conductance 
in units of e 2 //i, where e is the electronic change and h 
is Planck’s constant. This has relative uncertainties typ¬ 
ically smaller than 10 _1 ° between different samples and 
plateaus, which is an unprecedented level of accuracy for 
semiconductor systems which usually have unavoidable 
sources of disorder [l] . The origin of the precision is now 
understood to be due to the Hall conductance being a 
topological quantity related to the Chern number jH 3;]- 
Understanding topological states of matter continues to 
gain importance, where there is currently an intense ef¬ 
fort to investigating topological insulators Q and apply¬ 
ing these concepts to quantum computing using topolog¬ 
ical error correction codes, which have the highest error 
thresholds to date @, !G]. 

Realizing the QHE, and related topological states of 
matter, in systems other than semiconductors has there¬ 
fore become an important pursuit in several fields. For 
Bose-Einstein condensates (BECs), there is a well-known 
equivalence between a magnetic field and rotation that 
allows for applying a vector potential to charge neutral 
atoms [7, 0. This equivalence has suggested that the 
QHE may be accessed using BECs, when the bosons oc¬ 
cupy the lowest Landau level. Due to the interacting na¬ 
ture of the atoms in the BEC, it has also been predicted 
that the bosonic version of the fractional QHE should be 
observable, complete with non-Abelian quantum states. 
Experimentally, among technical challenges such as heat- 



FIG. 1: (Color online) Schematic configuration considered 
in this paper. (a)(b) A Bose-Einstein condensate (BEC) is 
confined in a long channel along the ^-direction with one or 
several vortices (each marked with a circle, rotation orienta¬ 
tion as marked) present in the central region. Depending on 
the position of the vortices, the condensate flows with a total 
current J x . Contour integrals for various paths C as discussed 
in the main text are marked. Dashed lines are contour in¬ 
tegrals that are at x = ±00 which contribute zero, (c) The 
experimental configuration considered in this paper. We as¬ 
sume that the density profile along the channel is independent 
of x, with the exception of density dips corresponding to vor¬ 
tex cores. Vortices are produced by optically based techniques 
such as stirring or giving angular momentum to the BEC with 
Laguerre-Gauss modes. Blue-detuned lasers pin the vortices 
and move them to various locations of the BEC. 


ing, the difficulties are to precisely match the rotation 
frequency to the trapping frequency, which has so far 
prevented realizing the QHE in cold atoms. In addition, 
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there has been a large amount of interest in realizing syn¬ 
thetic magnetic fields, and more generally, gauge fields 
[§]. Experimentally this was realized in Bose-Einstein 
condensates with the production of vortices without ro¬ 
tation 0- Several proposals for using such artificial 
magnetic fields to realize the fractional QHE 00, 
anomalous QHE 0], and quantum spin Hall effect 141 . 
Some of the remarkable progress towards realizing such 
schemes experimentally include the demonstration of the 
superfluid Hall effect 0], spin Hall effect 0, and mea¬ 
surement of the Chern number in Hofstadter bands ra¬ 
in this paper we present an alternative and very simple 
approach to observe integer QHE behavior in a BEC. Our 
scheme does not possess an strict mathematical equiva¬ 
lence to the QHE as the approach as described in Refs. 
[2, 8} (or the bosonic version of it). Nevertheless, it pos¬ 
sesses several essential characteristics in common. The 
quantization occurs with respect to the same observ¬ 
able as the standard QHE — the current response of 
the condensate — and occurs in units of m/2h, where 
to is the atomic mass. Furthermore, this quantization 
can be shown to originate from a topological quantity 
related to the phase of the wavefunction in the presence 
of vortices. In the standard QHE it is known that the 
quantization occurs due to presence of vortices in the 
Brillioun zone {2J, in this respect the origin of our effect 
is the same, except that the vortices are in real space. If 
we consider the origin of the vortices as rotation, we ob¬ 
taining a striking similarity to the standard conductance- 
magnetization quantization curve seen in the QHE, with 
discrete plateaus in the current response crossing over to 
linear behavior in the limit of many vortices. 


channel, it has no density dependence in the x-direction, 
i.e. 


ip(x,y) = V P(yy S(x ' y \ (!) 


where p(y) is a real function representing the density 
of the order parameter, and S(x,y) is the phase. No 
assumptions are made about the phase distribution. 
Strictly speaking, the presence of vortices already vio¬ 
lates the assumption of x- independence in the density as 
this implies local zeros in the condensate density. How¬ 
ever, if the area occupied by the vortices is small relative 
to the total area, we show in Appendix [A] that this is a 
negligible contribution to the total current. 

Using the fact that the current is j = ■^p'VS, the total 
current J x can then be evaluated to be 


JX 


h 

m 




dx 


dS(x, y) 
dx 


( 2 ) 


The crucial observation is that the integral in the x di¬ 
rection may be evaluated exactly, as the density is inde¬ 
pendent of x. This follows from the well-known general 
topological property of the velocity v(r ) = ^ VS(r ) in a 
BEC. For an arbitrary path around N vortices, the line 
integral is 0. 

Ivdl = — J)VS-dl = 2TrN—. (3) 

J m J m 


Now consider our particular geometry with a single vor¬ 
tex. For an arbitrary contour that extends to x = ±00 
such as that shown in Fig. [TJt, the contour integral can 
be written 


II. VORTEX DISPLACEMENT-CURRENT 
RELATION 

A. Single vortex configuration 

Figure Q] shows the basic configuration that we consider 
in this paper. The BEC is assumed to exist in a long 
channel, with the channel running in the x-direction. A 
strong ^-confinement allows us to consider an effectively 
two-dimensional system, such that the dynamics are en¬ 
tirely in the x-y plane. The a:-confinement should be very 
weak, such that a net current can freely flow in the re¬ 
direction. We furthermore assume that a small number 
of vortices are present in the center of the BEC. Let us 
now consider the total current flowing in the x-direction, 
equal to total velocity of all the atoms in the condensate 
in the x-direction, 

Jx = = f drj x (r) 

to J 

where m is the mass of the atoms, j(r) = 
['*/’( T ’)*' ( 7'0( r ) — V’( r )'^V , ( T ’)*], and V’( r ) is the order- 
parameter of the BEC. The primary assumption that we 
make is that due to the BEC being present in the long 


I Cb - I Ca = 27t (4) 

where 

I c = jvS-dl , (5) 

and C is a contour that runs from x € (— 00 , 00 ) but 
may take any path along the way. Ca is a contour that 
runs above the vortex, Cb runs below the vortex. The 
contributions from the edges at x = ±00 (dashed lines in 
Fig. [l}i) give zero contribution as these are perpendicular 
to the center of mass momentum kg of the BEC in the 
x-direction, and are very far from the vortex. From the 
symmetry of the configuration we deduce that 


Ic A = ko-n 

Ic B =k 0 + n, (6) 


where fco is the center of mass momentum of the conden¬ 
sate in the x-direction. For the specific case of a contour 
that is a straight line parallel to the x-axis, we have 



dx 


dS(x, y) 
dx 


= k 0 ± 7r, 


( 7 ) 
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where the sign is positive if y < y± and negative if y > y \, 
and yi is the y-coordinate of the vortex. 

We may now evaluate 0 directly to obtain 


j X 


HkoJ\f h 

—— + x— A y 
m 2m 


( 8 ) 


Here A f is the total number of particles in the BEC and 
the spatial asymmetry parameter is 


A y = dyp(y) - 


dyp{y)- 


(9) 


'y i 


The first term in Eq. ® is a trivial overall offset to 
the current and is independent of the vortex component. 
The second term however results entirely from the topo¬ 
logical phase of the vortex. The physics described by ® 
is rather simple: for various vortex displacements in the 
y-direction in a condensate, a current in the ^-direction 
proportional to the parameter ® is produced. A y is a 
parameter that counts the difference between the number 
of particles above and below the vortex. 


B. Multi-vortex configuration 


This argument is easily generalized to the case with 
multiple vortices. Using the labeling for the contour in¬ 
tegrations as shown in Fig. [l]r, we obtain a system of 
equations satisfying 

I Ck — I C i = 2(fe — l)n, (10) 

where k,l £ [0,1V], N is the number of vortices, and 
we have again assumed that all the vortices are located 
far away from the boundaries such that the contours at 
x = ±oo do not contribute. This can be solved to give 

I Cl =k 0 + (2 l - N)n. (11) 


This allows us to write the current-asymmetry relation 
more generally for the multi-vortex case 

J x = ± G xy Ay , (12) 

where we have defined the Hall conductance-like quantity 


G 


xy 



(13) 


It is clear that this is quantized in units of ^ for each 
vortex that is present. Here the asymmetry parameter is 


N 


A , = Y. 


21- n r m+i 


1=0 


N 


dyp(y) 


(14) 


where yu is the y-coordinate of the fcth vortex and we 
have defined yo = oo and yN+i = —oo. The power of a 
relation such as m - as is also true for the QHE in 
semiconductors — is that all the measurable quantities 
are easily accessible yet lead to a non-trivial quantum 
property of the system. 


C. Connections to Laughlin’s gauge argument 


In the previous sections we have derived a connection 
between the net current flowing in a BEC with displace¬ 
ments of vortices. This was derived from the topologi¬ 
cal integral of the phase around a vortex. This can be 
viewed also from the point of view of the Magnus force 
when moving the vortex 0, [13. For the homogenous 
case one may use an adaptation of Laughlin’s gauge ar¬ 
gument [2lT - [23l | to derive @ in a limiting case. 

Let us first derive the current-asymmetry relation for 
infinitesimal displacements of the vortex. Starting from 
® , consider moving the vortex from y 1 to y\ + 6y. The 
change in the current is 

d J x = ^ [A y {y 1 ± Sy) - A y {y 1 )] 

u ryi+Sy 

= — dyp(y) 

mj yi 

~ —dyp(yi)- (15) 

m 

Associating the local density n = p(y\) and the momen¬ 
tum J x = p x /m , we obtain the relation 

5p x ~ ZithnSy. (16) 

Using a similar argument one may derive the current for 
a multi-vortex configuration 

h 1 N 

SJ X « — Sy— p(yi )■ (17) 

m N 

1=1 


For the homogenous case with p(yi) = n this reduces to 

(USD- 

On the other hand, we may obtain a similar relation 
using a modification of Laughlin ’s g auge argument (see 
for example Sec. IXB of Ref. [2l| for a discussion of 
this). Extend the two dimensional x-y plane to a torus 
and consider threading a flux in the y direction. The di¬ 
mensions in the x and y directions are considered to be 
L x and L y respectively. This may be achieved by creat¬ 
ing a vortex-antivortex pair and moving them apart in 
the y direction until they annihilate at the opposite side 
of the torus. The movement of the vortices in the y di¬ 
rection eventually give a momentum to the whole system 
in the x direction, which is the effect we are interested 
in. The phase in the x direction that results from the 
flux threading is constrained to be e 27r “/ ix , i.e. the mo¬ 
mentum given to the system is 


A p x 


2 nhAf 


L 


X 


= 2irhnL y 


(18) 


where n = Yl • Integrating csD gives USD- 

While the Laughlin gauge argument gives qualitative 
agreement to our results, there are also several differ¬ 
ences. Firstly in Laughlin’s argument one generally con¬ 
siders the momentum difference before and after the flux 
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FIG. 2: (Color online) Stationary solutions of the Gross- 
Pitaevskii equation with a single vortex at r\ = (0, Sy ) with 
(a)(b)(c) 5y = 0 (d)(e)(f) 6y = 20. Plots for the conden¬ 
sate (a)(d) density, (b)(e) phase, (c)(f) momentum distribu¬ 
tion are shown. Length scales are in units of the healing 

length £ = 2 mgn where n is the maximum density of the 
BEC, along y = 0. We assume parameters tkj y /Eo = 0.02, 
Vo/Eo = 1, fco = 0, where the energy scale is Eo = 2 '' n ^ = gn. 
The delta potential is located at the vortex position r i. 


is threaded. Our relations describe the current relation 
for an arbitrary vortex configuration. This is desirable 
particularly when a strict proportionality would like to 
be extracted, as is our case where m/h is the quantity 
that requires estimation. Secondly, our results do not 
assume homogeneity in the y-direction, which may be 
important to include in realistic BECs. 


III. NUMERICAL SIMULATIONS 


To illustrate the effect, we perform numerical simula¬ 
tions of the Gross-Pitaevskii equation 


,/lv 

,h Tt = 


fi 2 V 2 1 


2m 


muly 2 


Vo^6(r - r k ) + g\ip\ 2 ip, 
k 

(19) 


where oj v is the trapping frequency in the y-direction, 
Vo is the strength of the pinning potentials at locations 


rk, and g is the interaction strength. Using a real space 
second order finite difference time decimation (FDTD) 
method on a 100 x 100 site grid simulating the region 
x S [—x max , x max ], y £ [-y max , y max ] evolving in time using 
the backwards Euler method. In order to obtain station¬ 
ary vortex states we evolve in imaginary time starting 
from an approximate vortex wavefunction. In order to 
have stable vortex solutions under imaginary time evo¬ 
lution, we pin the vortex using a local delta function at 
the desired vortex position, which has a negligible effect 
on the phase and density of the condensate wavefunc¬ 
tion. To support a constant current in the x-direction, 
we employ (anti-)periodic Mobius boundary conditions 
ip(-Xm**,y) = (±1 ) JV '0(x max , —y) for an even (odd) num¬ 
ber of vortices. The phase factor of (±l) w is necessary 
as each vortex flips the phase by either 7r going from 
x = — oo -A oo, as discussed in the main text. The 
Mobius boundary conditions are necessary as the cur¬ 
rent far to the left of the vortices are predominantly in 
the +y direction, whereas far to the right the current is 
in the —y direction, or vice versa. Parameters are chosen 
consistent with BECs with large particle numbers where 
the condensate radius is much larger than the healing 
length gjj . 

Our results for a single vortex are shown in Fig. [2j The 
real space, phase, and momentum space distributions for 
two vortex positions are shown. In both cases the imag¬ 
inary time evolution ensures a quiescent stationary state 
of the BEC which smooths out the density fluctuations 
in the x-direction. The lack of such density variations is 
our primary assumption, and is the corresponding situa¬ 
tion in the standard QHE to an equilibrium state free of 
transient dynamics. The real space images [Fig. [2fa)(d)] 
show that such a state is achieved with a stable vortex 
pinned at their respective positions. The phase variations 
[Fig. [2][b)(e)] agree well with the discussion relating to 
Fig. [lja), showing that there is a phase change of ±7r 
depending on whether the y-position is above or below 
the vortex. For the central vortex position the momen¬ 
tum distribution is symmetrically distributed such that 
the average momentum (and hence current) is ( k x ) ~ 0. 
For the displaced vortex position, the momentum distri¬ 
bution shifts to the right, indicating a non-zero ( k x ) > 0, 
as predicted by the relation ©. We emphasize here that 
the center-of-mass momentum is k$ = 0, so that the net 
momentum in the x-direction in dSJ) results entirely from 
the vortex displacement. 

Similar results are obtained for the multi-vortex case as 
shown in Fig. [3] The same procedure as the single vortex 
case is repeated for a configuration of three vortices in an 
equilateral triangle, displaced by various y positions. The 
phase relation m can be seen to hold by taking lines at 
various y-positions, giving phase shifts of —37 t, —7t, 7r, 37r. 
The momentum distributions again shift towards the pos¬ 
itive ( k x ) > 0 direction once the vortices are displaced. 
The current-asymmetry relation for various numbers of 
vortices and y-displacements are shown in Fig. QJi. We 
see a perfect proportionality relation as predicted by m- 
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FIG. 3: (Color online) Stationary solutions of the Gross- 
Pitaevskii equation with three vortices with coordinates ri = 
(1(V3,20 + Sy), r 2 = (-10^3 ,Sy), r 3 = (10^-20 + Sy), 
where (a)(b)(d) 5y = 20 and (c) 8y = 0. Plots for the conden¬ 
sate (a) density, (b) phase, and (c)(d) momentum distribution 
are shown. The same parameters and units as Fig. [2]are used. 


in the rotating frame. Using this we may then calcu¬ 
late a conductance quantity G xy as the vortex positions 
are changed. Taking this value as the vortex number in 
G xy , we obtain a curve which is remarkably reminiscent 
of the conductance-magnetic field relation in the QHE 
(Fig. SKb)). For slow rotations we recover the conduc¬ 
tance plateaus corresponding to low vortex numbers. For 
fast rotations there is a proliferation of vortices, and si¬ 
multaneously the resistance R xy diminishes as oc 1/JV, 
which gives a linear relation. We note that a similar 
curve can also be obtained by plotting 1/fl versus R xy , 
although this gives a square root relation as the critical 
Q/u} c —> 1 is approached. In practice, it is likely that 
directly rotating the BEC is not the best experimental 
method for vortex generation (see Sec. IIVI) . Neverthe¬ 
less, as rotations are the corresponding quantity to the 
magnetic field in the BEC case, the equivalence to the 
standard QHE is most clearly illustrated by the conduc¬ 
tance quantization as shown in Fig. |U(b). 


IV. EXPERIMENTAL REALIZATION 


A linear fit to the data points give quantization to the 
integer multiples of Go = h/2m. While good agreement 
with the theory is observed, we attribute discrepancies to 
exact integral quantization to the relatively short chan¬ 
nel that we use in the simulations, of which the length is 
only of the order of the width. The very precise quanti¬ 
zation as seen in the QHE is a result of the whole system 
contributing to the conductance. Thus as the system size 
is increased we expect the precision to improve (see also 
Appendix EJ. In experimental systems we expect that 
much longer channels can be produced, and so can ben¬ 
efit from this scaling effect. As our arguments are based 
on rather general topological considerations of the phase 
around vortices in a BEC, the effects of local disorder are 
not detrimental to the effect, assuming that they do not 
cause large variations in the density along the channel. 

We now show the explicit quantization relation of the 
inverse conductance R xy = 1/G xy to show most clearly 
the analogy with the QHE. Let us assume that the vor¬ 
tices are originally produced by a rotation of frequency 
Cl. The number of vortices that are produced in two di¬ 
mensions can then be estimated to be [25j 


N = 


Q/uJc 


( 20 ) 


where w c is the critical frequency which gives rise to 
the proliferation of vortices, n is a dimensionless pro¬ 
portionality constant, and the bracket [] rounds to the 
nearest integer. Let us now consider that in this situa¬ 
tion we examine the J x current and vortex positions yi 


We finally discuss the likely experimental configura¬ 
tion of our proposal. A oblate atomic BEC with trap¬ 
ping frequencies satisfying w 2 > to y uj x would be pre¬ 
pared, such as the dynamics would primarily occur in 
the x-y directions. To ensure that such vortices can ex¬ 
ist in the channel, the y-confinement and density should 
be such that the width of the BEC is larger than the 
healing length. The schematic configuration is shown 
in Fig. me). Starting from such a configuration, small 
numbers of vortices would be generated. This is most 
suitably done with optically based techniques such as 
stirring with a blue-detuned laser [2(il - [28| . or adiabati- 
cally introducing angular momentum to the BEC using 
Laguerre-Gauss modes f29| - f3T1 [ . In order to extract G xy 
from the current and asymmetry relation, a variety of 
vortex positions are required. To achieve this in the most 
controlled fashion, pinning of the vortices at a given loca¬ 
tion is desirable. Vortex-antivortex pairs may be reliably 
generated and pinned at a desired location by the use of 
two blue-detuned lasers moved through the BEC [H, . 

Once the vortices are generated, the density distribu¬ 
tion of the BEC is estimated using high-resolution spatial 
imaging techniques such as scanning electron microscopy 
P . Other in-situ methods such as as phase contrast 
imaging could be used to obtain a density distribution of 
the BEC to identify the vortex positions [35, 36|. The 
total current J x = -^ k k x \ipka,k v | 2 can be extracted 
from the high-resolution velocity distribution of the BEC, 
which can performed by Bragg spectroscopy [37J or time- 
of-flight imaging [3S| . 

In order to measure the conductance plateaus, one may 
directly use the relation m, or alternatively the differ- 
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FIG. 4: (Color online) (a) Total current J x versus asymmetry 
parameter A y for various vortex number N. Points show the 
numerically evaluated values and lines are fits to the data. 
The current J x is calculated in units of Go = the ele¬ 
mentary unit of conductance, such that the gradient should 
be an integer. A linear fit to the data gives a gradient of 
0.99, 2.00, 2.99 for the N = 1,2,3 vortex cases respectively, 
(b) Resistance R xy = 1/G xy versus the rotation parameter 
\/l — (fi/u>c) 2 - Resistance is measured in units of inverse Go 
and Ac = 1 is used. The number of vortices generated by the 
rotation are labeled. Plateaus are rounded to account for the 
uncertainty in vortex number between plateaus. 


ential form CSD which can be written as 
h f yi+5y 

SJ X = — / dyp(y) 

mjyi 

~ —Syp(yi) ( 21 ) 

m 

where the current per particle is 


SJ X 


<4(2/1 + Sy) - 4(2/i) 

a r 


and the normalized density is 


p(y) 


p{y) 

M ' 


( 22 ) 


(23) 


The advantage of writing the relations using normalized 
quantities is that they become insensitive to number fluc¬ 
tuations in the BEC which may originate from variations 
in the initial conditions or particle loss. Thus during the 
measurement process, only the relative position and mo¬ 
mentum distributions are required. The ratio h/m is 
then extracted by looking at the differential variation of 
the current 5J X with variations of the vortex position Sy. 
The spatial measurement then plays the role of obtain¬ 
ing p(yi), i.e. the relative density of the condensate at 
2 /-coordinate of the vortex. This can also be obtained by 
interpolating the density without the vortex core. 


able to measure the absolute mass of the atoms, rather 
than the relative atomic mass. Currently Penning traps 
are the most precise mass spectrometers [3E| achieving 
relative uncertainties of typically ~ 10 -1 ° for the relative 
atom mass. The (absolute) atomic mass unit is known 
to an larger uncertainty of ~ 10 -8 jdcj. If an extremely 
precise measurement of absolute mass of an atom became 
possible, this would allow for a redefinition of the kilo¬ 
gram as a fixed number of atoms of a particular type; 
87 Rb for instance (see Appendix 0. Since the conduc¬ 
tance Go is the quantity that would be measured in our 
proposal, the mass would be measured in units of the 
Planck’s constant. This is consistent with other meth¬ 
ods that aim to contribute to a redefinition of the kilo¬ 
gram such as the watt balance and silicon atom counting 
methods. As another potential application, the quantity 
h/m itself is of interest in the context of determination 
of the fine structure constant, by combining estimates of 
the Rydberg constant, the relative mass of an atom, and 
h/m [401 - The first two factors are determined to better 
than 10 -10 , while currently h/m can be estimated to a 
level of 10 -8 . 

The spatial and velocity distribution measurement 
methods both have a finite resolution which may appear 
to severely limit the precision of the quantities to be esti¬ 
mated in (ET1) . However, bulk quantities such as 8J X scale 
well with finite resolution, e.g. for Simpson’s rule the er¬ 
rors scale as the fourth power of the discretization. For 
p{yi), assuming that the vortex is pinned by the blue- 
detuned laser so that it is always present at the same 
location, then one can measure the relative density with 
high accuracy by many repetitions of the experiment. 
To minimize the uncertainty of the density measurement 
it is advantageous to move the vortex in the vicinity of 
the maximum density of the BEC, where the derivative 
with position is zero, i.e. p(y\ ± Sy) ~ p{yi) ± Sy. In 
this way errors due to finite resolution can be mitigated. 
This is also advantageous in terms of corrections to the 
current due to the vortex (see Appendix lAl) . Finally, Sy 
can be set by the pinning laser which is controlled and 
hence is not measured directly in the BEC. We therefore 
estimate that the main sources of error will arise from 
density and thermal fluctuations in the BEC, and other 
experimental issues such as calibration of measured and 
controlled quantities, and repeatability. Although it is ul¬ 
timately an experimental question of whether sufficiently 
low uncertainties can be attained, at the level of the sys¬ 
tem the topological nature of the observable should make 
the measurement of G xy rather robust under a variety of 
conditions, in analogy with the standard QHE. 


V. APPLICATIONS 

Since the conductance G xy is quantized in units of Go, 
our result presents the possibility of an novel method of 
measuring the mass of atoms in the condensate. We note 
that this would be a mass spectrometer that would be 


VI. SUMMARY AND CONCLUSIONS 

In summary, we have shown an alternative method of 
investigating integer quantum Hall physics using BECs 
with small numbers of vortices. While this approach does 
not have a precise mathematical equivalence of previous 
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approaches 00 , much of the essential physics is in com¬ 
mon, where the topological phase around the vortices 
give rise to a quantized Hall conductance-like quantity. 
In both cases the current response along the channel is 
measured, but here the proportionality is with respect 
to an asymmetry parameter, as opposed to the standard 
QHE where it is the potential difference in the trans¬ 
verse direction. Alternatively the density at the same y- 
coordinate as the vortex can be measured instead of the 
asymmetry parameter. A potential application is to use 
the quantization of the conductance in units of h/2m as 
a novel way of measuring the absolute mass of the atoms, 
using a topological invariant of the system. While in this 
paper we have implicitly assumed an atomic BEC, other 
types of BEC, such as exciton-polaritons [UHU, should 
also be suitable to observe the effect as the effect relies 
only on the topological phase of the condensate. 
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Appendix A: Corrections to the conductance 
quantization due to vortices 

The main assumption made in our calculations is Eq. 
(3), that the BEC’s density is uniform in ^-direction. 
The presence of vortices clearly violates this assumption. 
In this section we show to what extent the presence of 
vortices make to the final result. 

The condensate wavefunction with a single vortex can 
be written 


V’O, V) = [f(y) ~ A /(r v )] e lS(x ’ y \ (Al) 

where A/(r v ) is the density deviation due to the presence 
of the vortex, f(y) = p(y), r v = |r —ri|, and the vortex 

is located at the position r±. This is typically a positive 
quantity that has a maximum at the vortex core and 
approaches zero at a distance of the order of the healing 
length. The current can be written 

Jx =— f dxdy ^ [f 2 (y) - 2/(y)A/(r v ) + (A/(r v )) 2 ] . 
m J ox 

(A2) 


Evaluating this expression as in the main text, we obtain 


JX 


h 

m 


koAf J— A + A J x 


(A3) 


where we have used the same definitions as the main text 
and 

r ae 

A J x = J dxdy—Af(r v ) [A/(r v ) - 2 f{y)} . (A4) 

The main contribution to the integral is the region 
around the vortex, which is assumed to be a small region 
in comparison to the scale of the function f(y). Shifting 
the coordinates to the center of the vortex, we can then 
approximate the phase in this region as being 
(here Sq is a constant), and f(y) ~ /o + fiy + ■ ■ ■ to first 
order. This gives = —2^, and hence 


A J x = — J dr v (A/(ry)) 2 j d^sin</> 

+ 2/ 0 J dr v A/(r v ) J d(j) sin</> 

+ 2/i J dr- v r v A/(r v ) J d</>sin 2 (j) + ... (A5) 


The first and second terms are zero due to the integral 
over the phase. Thus for local densities f(y) that are 
flat give zero correction to the current. The first order 
correction enters when there is a gradient /i in the local 
density. The last term can be evaluated and we obtain 

A J x = 27 t/i J dr v r v A/(r v ) + ... (A6) 

Let us now estimate the order of magnitude of each of 
the terms in (IA3D . Let us write the average density of 
the BEC as n ~ > w h ere L x ,L y are the lengths of 

the BEC in the x , y directions. Taking the magnitudes 
of /o ~ y/n 7 /i ~ y/n/l, and A/(r v ) ~ -y/n, where l is the 
length scale associated with the gradient, we have 


hn 

m 


koL X Ly + Ly + 


A. 


(A7) 


where A vorte x is the area of the vortex. For the correction 
due to current due to the vortex to be negligible, we thus 
require that 

L y » ^1. (A8) 

This is satisfied if the area of the vortex is very small 
compared to the condensate, and if the vortex is present 
in a very flat region of the BEC. 


Appendix B: Mass standards 

The aim of a mass standard is to create a new defi¬ 
nition of the kilogram, which is currently defined as the 
mass of a platinum-iridium alloy artifact stored in France 
iTT*H . As with other definitions such as the meter and the 
second, it is desirable to use fundamental constants of 
nature rather than material artifacts, for several reasons 










such as stability, reproducibility across the world, and 
other practical issues causing the artifact’s mass to drift 
in time. We give some more detail on how our proposal 
would be connected with the kilogram mass standard. 

One approach to a mass standard is to define 1kg to be 
equal to a certain number of carbon-12 atoms. However, 
there is already the definition that 1 mole (= Na, Avo- 
gadro’s number) of carbon-12 is 12 grams exactly. This 
means that when the kilogram is redefined the constraint 

N A m c = 0.012 (Bl) 

must also be satisfied, where me is the mass of 1 carbon- 
12 atom (in kg). Currently, both N A and me are ex¬ 
perimentally determined quantities. In a redefinition of 
the kilogram using carbon atoms, the Avogadro con¬ 
stant would be fixed to a particular number, for example 
N a = 6.02 2 1 4 x 10 23 exactly. Then according to the 
constraint (ED this fixes me also exactly. Alternatively, 


me could also be fixed, and this would fix N A according 
to ED- Fixing me or Na is equivalent in this sense. 

A very precise measurement of the mass of the atoms 
can therefore equivalently contribute towards the mass 
standard. While a BEC of carbon-12 is not practical, 
the relation ED can be converted to something more 
convenient for this purpose. Introducing the mass of the 
atom species that undergoes BEC m (such as 87 Rb) we 
have 

N A —m = 0.012. (B2) 

m 

The relative mass ratio can be measured very pre¬ 
cisely (typically 10~ 10 relative error) using Penning traps 
fiH. In this way by measuring the (absolute) mass of the 
BEC atom, this can be used to determine the Avogradro 
constant, which defines the mass standard. 
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